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ABSTRACT 

We investigate how the dynamical state of a turbulently supported, 1000 M Q , 
molecular cloud affects the properties of the cluster it forms, focusing our discussion 
on the star formation efficiency (SFE) and the initial mass function (IMF) . A variety 
of initial energy states are examined in this paper, ranging from clouds with £g ra v| = 
0.1-Ekin to clouds with \E. 
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1 INTRODUCTION 



The internal motions in giant molecular clouds (GM Cs) are 
dominated by supersonic turbulence (for example, Larson! 
19811; iMversI Il983l; iDame et al.l Il986l ; ISolomon etall 1 19871 ; 



Hever fc Bruntl 2004 ). These m otions leave their imprint 



on the cloud's internal str ucture (|Mac Low fc Klessen|[200l ; 
lElmegreen fc Scald |2004| ). resulting in dense regions, of- 
ten referred to as clumps or cores, that are observed to 
have a wide range o f dynamical states (|Ballv et al I I 19871 ; 
IWilliams et al.|[l994l ). The aim of this paper is to explore 
how the initial dynamical state of star forming regions can 
influence both their final stellar mass functions and their 
star formation efficiency (SFE). 

IClark fc Bonnelll <|2004l ) suggested that stars can form 
in dynamically unbound molecular regions, where the in- 
ternal turbulent motions are in excess of the region's grav- 
itational energy. Depending on the scale, such regions can 
be star forming cores or clouds, and could thus give rise 
to a wide range of cluster sizes. For dynamically unbound 
regions to be capable of producing stars, they need to 
contain multiple thermal Jeans masses, or undergo bet- 
ter than i sothermal cooling during the shock compression 
phase (e.g. lChapman et al.lll992l ; lBhattal et al 1ll998l ; lLarsonl 
120051 ). This is simply because isothermal shock compression 
(in 1-D) does not significantly increase the number of Jeans 



10 -Ekin j and for both isothermal and piece- wise 
polytropic equations of state (similar to that suggested by Larson). It is found that 
arbitrary star formation efficiencies are possible, with strongly unbound clouds yielding 
very low star formation efficiencies. We suggest that the low star formation efficiency 
in the Maddelena cloud may be a consequence of the relatively unbound state of its 
internal structure. It is also found that competitive accretion results in the observed 
IMF when the clouds have initial energy states of |i? gra v| ^ £-kin- We show that 
under such conditions the shape of the IMF is independent of time in the calculations. 
This demonstrates that the global accretion process can be terminated at any stage 
in the cluster's evolution, while still yielding a distribution of stellar masses that is 
consistent with the observed IMF. As the clouds become progressively more unbound, 
competitive accretion is less important and the protostellar mass function flattens. 
These results predict that molecular clouds should be permeated with a distributed 
population of stars that follow a flatter than Salpeter IMF. 

masses in the shock layer dElmegreen fc Elmegreenl [l978l ; 
|Poroshkevicn]|l980l ; iLubow fc Pringld Il993l ; IClarkej ll999l ). 
This fact is apparent in the numerical studies of turbulent 
molecular clouds, which show that the number of stars is 
always comparable to the number of initial Jeans masses in 
the simulation, despite the huge densit y contrast produced 
by the turbulence or colliding flows (IKlessen et al 1 12000| ; 
Icllttins et al.ll2003l ; IClark fc Bonnel]||2005l ; lLarsonll2007r i. 



Dynamically unbound star forming regions have 
three important features. First, star formation occurs 
on a similar timescale to the dispersion of the region. 
Such a pictu re is consistent with the rapid star forma- 



tion model |Vazquez-Semadeni. Passot fc Pouquej 



Ballcstcros-Paredes. Hartmann fc Vazq ucz- Scmadcni 



199" 



19!)!; 



Elmegreenl |2000|; lHartmann. Ballesteros-Paredes fc Berginl 



2001 



Pringle. Allen fc Lubo^l l200ll ). This 



pro- 

cess has gained r ecent support from the studies of 
iGlover fc Mac Low! (|2007al lbh. which show that H 2 forma- 
tion can occur on a dynamical time in both turbulent or 
freely collapsing gas. The second feature is that unbound 
clouds can naturally result in lower star formation efficien- 
cies than their bound counterparts, simply because much 
of the gas is unable to rid itself of the excess kinetic energy 
and achieve a bound state. The efficiencie s can range from 
50 percent at protostellar core scales (|Clark fc Bonnelll 
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Figure 1. Evolution of the column density distribution in the gas for 2 of the a k i n = Ekin / |S grav | = 2 simulations, for an isothermal 
(top) and a piece- wise barotropic EOS (bottom), similar to that discussed by Larson (2005). Both the simulations have the same initial 
turbulent field. The colour table in the images is stretched (logarithmically) over 5 orders of magnitude in column density, ranging from 
0.001 to 100 g cm -2 . 



|2004| ). to b etween 5 and 10 percent at the larger scales 
of a GMC (jClark et alj|2005l ). Finally, IClark et al.1 (|2005l ) 
also found that at the GMC scale, the unbound cloud can 
produce a series of clusters that reproduce some of the 
features found in OB sub-groups. 

The aim of this present study is two-fold. First, 
we wish to explore how the initial energy balance 
in the cloud can affect the star format io n effi ciency, 
extending t he w ork of IClark fc Bonnelll (12004 ) and 
Clark et al.l (120051) (see also Klessen. Heitsch fc Mac Low! 



200d. 



Heitsch. Mac Low fc Klessen 



Vazquez-Semadeni. Ballesteros-Paredes fc Klessen! 



Secondly, we wish to look at how the mass function of the 
protostars is influenc ed by such an environment. When 
competitive accretion (|Zinneckerll 19821 ; iBonnell et aLll2001al ; 
iBonnell et ail l2001bh is the dominant process controlling 
the distribution of stellar masses, it has been shown that 
the break in the protostellar mass function, and thus the 
characteristic mass , depends on the mean Jeans mass 
dKlessen et al.lll998l;lKlessen fc Burkertll2000l. I200H; iKlessenl 
l200ll ; Ijappsen et all 120051 ; IBonnell. Clarke fc Bate! l200ef )7 
For clouds which are out of global dynamic stability, the 
mean density (and thus the Jeans mass) are constantly 
changing. As such, the ability of competitive accretion to 
produce an IMF-like protostellar mass function under these 
conditions is not clear. 

To address these questions we perform a series 
of numerical simulations which follow the evolution of 
1000 M Q clouds, in which the initial ratio of turbulent ki- 
netic energy to gravitational potential energy is systemati- 
cal ly varied. These clouds , which are similar to those studied 
by IBonnell. Bate fc Vine] (|2005t ), can be thought of as over- 
densities in a larger GMC, which still contain residual kinetic 



energy from the supersonic flows that created them. Further 
details of the calculations are given in Section [2] and we dis- 
cuss the star formation efficiency and mass functions from 
these clouds in sections [3] and [4] respectively. The implica- 
tions of our results for star formation in GMCs in general are 
presented in Section [5] and a brief summary of our findings 
is given in Section [6] 



2 THE SIMULATIONS 

2.1 Computational Method and Cloud set-up 

The evolution of the gas was modelled using the Smoothed 
Particle Hydrodynamics (SPH) me thod (for a description 
of this technique, see for example iBenj Il990l . iMonaghanl 
1 1991 lMonagharJl2005l ) , and the simulations we present here 
were performed on the UK Astrophysical Fluids Facility 
(UKAFF), a collection of IBM Power 5 nodes. The code 
utilises adaptive smoothing lengths for the gas particles, al- 
lowing high density regions to be followed with increased 
resolution. Gravity is calculated using a binary tree algo- 
rithm, which also serves to construct SPH neighbour lists. 

The clouds in this study have a mass of 1000 Mq , which 
we model using 2 x 10 6 SPH particles in an initially uni- 
form density sphere (see below for the initial density values). 
With this set-up, the simulations have a mass resolution 
of 0.05 Mq, sufficient to capture the frag mentation of the 
gas t o below the hydrogen burning limit (|Bate fc Burkert] 
1 1997ft . The clouds are also given a turbulent velocity field, 
which is allowed to decay freely via shocks. The turbulent 
velocity fields imposed here are described by a power spec- 
trum of P(k) oc fc -4 , which is consistent wi th the turbu- 
lent motions observed in molecular clouds (|Larsonl Il98ll ; 
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iMversfc Gambia [19991: iHever fc Bruntl [20o3). In this pa- 
per, we investigate the systematic changes that occur in the 
properties of the star formation as the dynamical state of 
a cloud is altered. We therefore use the same realisation of 
the turbulent field for most of the simulations. We do how- 
ever perform several simulations with a different turbulent 
field, to check that any systematic trends in the star for- 
mation properties are not due to any specific feature of the 
initial turbulent field. To adjust the dynamical state of the 
clouds, we change the initial ratio of kinetic to gravitational 
potential energy, akin = -Bkin / |S grav |, to vary from 10 to 
0.1. 

We adopt two types of equation of state (EOS) in this 
study: the first assumes the gas is isothermal and the sec- 
ond is a barotropic EOS, similar to that proposed by |Larso"nl 
|2005l ). In the isothermal gas, we set up the gas to have 
a Jeans mass of roughly 1 M© , such that the cloud has 
initially 1000 Jeans masses. In these isothermal runs, the 
clouds have a temperature of 10 K, and start at a density of 
1.45 x 10 19 g cm -3 . For the barotropic EOS, we use 



T — T c — , for p < p c 

Pc 

and then, 

T = T C , for p > p c , 



(1) 



(2) 



where p c = 5.5 x 10~ 19 g cm~ 3 and T c = 7.5 K. The values for 
T c and p c , ensure that at the point where the barotropic EOS 
enters the isothermal regime (that is, at densities higher 
than p c ), the means Jeans mass is roughly the same as the 
initial value (roughly 1 Mq ) found in the pure isothermal 
clouds. The barotropic EOS clouds start at a density of 
1.83 x 10~ 20 g cm~ 3 and a temperature of 17.6 K. 

To model the star fo rmation in this stu dy we use 'sink' 
particles, as described bv lBate et alj (|l995T ), which involves 
replacing dense, self-gravitating, regions of gas with a point 
mass. Sink particles are capable of accreting further material 
and they can interact with the other particles in the simu- 
lation via gravity. The sinks are formed in these simulations 
once a particle and it's neighbours are bound, collapsing, 
and within an accretion radius, h acCl which is taken here to 
be 412 AU. The sinks are then able to accrete any bound 
SPH particles which fall inside the accretion radius. The 
gravitational interactions between the sink particles and all 
other particles in the simulation are also smoothed to h acc . 
Due to the size of the sink particles in this study, it should 
be stressed that their mass functions presented in this paper 
should not be thought of as a stellar object IMF. Instead, we 
suggest t hat they are more representative of a 'system' mass 
function (|Kroupa. Gilmore fc Toutlll99ll ; lKroupall2002ft . As 
such, we will refer to them simply as 'sink' mass functions 
in the remainder of this paper. 



2.2 General Evolutionary Properties 

For the clouds with akin ^ 1, we run the simulations for 
two initial free-fall times. In our discussion of these simu- 
lations, we will use the properties at the end of the run. 
For our weakly supported clouds, which have akin = 0.1, we 
only follow the evolution for 1 initial free-fall time. In these 
latter simulations, the extreme dynamics and environment 
make them computationally very expensive. However, since 



they also evolve faster, one gains enough information for our 
present study. 

In Figure [T] we show column density snapshots from 
two of the simulations in this study. Both have the same 
turbulent velocity seed, and have akin = 2. The evolution of 
the isothermal cloud is shown in the upper panels, and the 
evolution of a cloud with a barotropic EOS is shown in the 
bottom panels. The expansion of the cloud is apparent in the 
images, and one can see that the clusters move apart as they 
form, in contrast to the merging process that has been doc- 
umented in studies where akin = 1 initially (jBonnell et al.l 
l2003l ; lBonnell. Vine fc Batell2004h . One can also see that the 
two different EOSs look very similar, since the sites of star 
formation are selected by the turbulent velocity field. 



3 THE STAR FORMATION EFFICIENCY 

We will first look at how the star formation efficiency (SFE) 
in a molecular region can be controlled by its dynamical 
state. We define the SFE in this study as, 



SFE: 



M sin k 



(Af gas + Ms inks ) 



(3) 



where M s i n ks and M gaa are the total mass in sink par- 
ticles and gas, respectively. In Figure [2] we plot the SFE 
as function of time for all the simulations with a^n = 
Skin / |-Egrav| ^ 1. The isothermal calculations are on the 
left and the barotropic EOS are on the right. We plot time 
in units of both the initial free-fall time of the clouds, and 
in years. The free-fall time is 1.75 x 10 yr in the isothermal 
runs, and 4.92 x 10 5 yr in the barotropic EOS. 

The plots show that a wide range of SFEs are possi- 
ble, depending on the initial dynamical state of the cloud, 
with the SFE decreasing with increasing akin- These cal- 
culations also suggest that star formation is able to take 
place under a variety of conditions, provided that the region 
has a sufficient number of Jeans masses at the outset. While 
the SFEs in our simulations are clearly not converged, and 
without feedback one would not expect them to be, they 
do demonstrate a clear trend with the initial level of tur- 
bulent support. One must also consider here that the SFEs 
given in Figure [2] are upper limits to the true value, since 
sink particles do not contain any model of the feedback that 
may be important in inhibiting further accretion. Also, not 
all the gas accreted by the sinks would be expected to find 
its way into stars, espec ially for the high-mass objects (see 
IZinnecker fc Yorkdl2007l for further discussions). 

The barotropic and isothermal EOSs give very similar 
results, in terms of the SFE evolution as a function of the 
initial free-fall time. This suggests that the regions for col- 
lapse are selected for fragmentation at similar times in the 
evolution, and that these regions proceed directly to col- 
lapse on the corresponding free-fall time. The similarity of 
both the cloud structure and cluster distribution between 
the two different EOS (Figure [T]) also suggests that the gen- 
eral morphology of the region is controlled primarily by the 
turbulent motions flows, via shock compression and dissipa- 
tion of kinetic energy. It is also interesting to note just how 
quickly the SFE drops as the value of Okm increases. Given 
that most th e structure in a GMC is typically thought to 
be unbound (jWilliams et al.lll994h . one would expect the 
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Figure 2. The star formation efficiency, SFE = M s i n k s / (M gas + M^nks), as a function of time for the simulations with akin = 
^kin / l^gravl 1? 1. Two different random realisations of the initial turbulent velocity field were performed for a^in =1, 2, Sz 3.3 in the 
isothermal case, and for a^in = 2 in the barotropic case. Time is given in both units of a free-fall in the bottom x-axis and in years in 
the top x-axis. 



majority of star formation in GMCs to only occur in a few 
regions, in which the local SFE can be quite high. 

The weakly supported clouds, naturally, have a much 
higher rate of star formation. We plot in Figure [3] the 
SFE as function of time (free-fall units) for the clouds with 
ttkin = 0.1, and for comparison, we also show the SFE curve 
from the isothermal akin = 1 calculation. Clearly, these 
weakly supported clouds evolve much faster than their more 
strongly supported counterparts, reaching 3-4 times the ef- 
ficiency in only half the time. 



4 THE SINK PARTICLE MASS FUNCTIONS 

As mentioned in Section [1] above, the effect of the compet- 
itive accretion process on the IMF in dynamically unbound 
simulations is still unclear, and needs to be examined. In 
Figure [4] we plot the sink particle mass function for all the 
simulations that have a^ n = Skin / l-Egravl ^ 1. The isother- 
mal clouds are shown on the left and the barotropic EOS 
clouds on the right. Some properties of the sink particles at 
the end of simulations are given in Table [3] 

The panels on the left-hand side show the sink mass 
functions for the calculations with akin = 1. These have a 
simila r initi a l set- up to those performed by iBonnell et al.l 
l|2003l . 120041 . 120061 ). where it was demonstrated that for a 
mean initial Jeans of roughly IMq, the resulting sink mass 
function is similar t o the field-star IMF (for the driven tur- 
bulent case, sec also Klcsscn 2001). Our simulations yield a 
similar result. Although the mean initial Jeans mass is much 
higher in the barotropic EOS simulation, the final result is 
very similar to that from the isothermal calculation. As has 
been discussed by o thers (|Larson|l2005l ; [jappsen et al.ll2005l ; 
IBonnell et al.l [20061 ) . this is simply because the Jeans mass 
at the critical density, p c , where fragmentation can occur, is 
similar to that in the isothermal case. 



One can see from Figure [4] that there is a distinct trend 
in the sink mass functions as the value of akin is increased: 
the mass functions become flatter as the clouds become pro- 
gressively more unbound. A similar behaviour is seen in cal- 
culations of dr iven turbulence with decreasing driving scale 
|KIessenll200~lh . This feature is also seen in the barotropic 
EOS calculations, suggesting that the trend in the mass 
function is not simply due to the changing Jeans mass in 
the unbound clouds. One can see this trend in the prop- 
erties of the sinks, given in table [4] The quantity /m> im 
denotes the fraction of the total mass in sink pa rticles that is 
made up of sinks with masses IMq or greater. iMever et al.l 
(|2000l ) performed a similar analysis on a number of local 
embedded clusters (using the number of stars, rather than 
the mass). One can see here that the fraction of the mass 
function comprising sinks of 1 Mq or more is clearly rising 
as akin is increased, indicating an overall flattening of the 
mass functon. 

There are two effects which will play an important role 
in shaping the sink mass function in the unbound clouds: 
the timescale of interactions, and the degree of fragmenta- 
tion. First, we consider the timescale on which newly formed 
objects interact with their siblings. In the a^ n = 1 simula- 
tions, the clouds are characterised by a particular density. 
For the isothermal clouds, this is the initial mean cloud den- 
sity, while in the barotropic EOS clouds, this is density, p c , at 
which the cooling stops. When fragmentation occurs at this 
density, the timescale for neighbouring fragments to interact 
is the same as the timescale for each fragment to collapse to 
form a sink particle. Naturally, this is just the free-fall time 
associated with the density where the fragmentation sets in. 
In contrast, the background density in the unbound simu- 
lations is constantly decreasing. The regions which become 
bound, and so manage to form sink particles, decouple from 
the surrounding material, which is evolving on progressively 
longer timescales. As a result, the star-forming regions which 



The SFE and its relation to variations in the IMF 5 



Table 1. Wc give here the properties of the sink particles after 
2 initial free fall times have elapsed, for all the simulations with 
Q kin ^ 1- For comparison, we also show the results from one of 
the simulations with a^ in = 0.1, with the results taken after 1 
initial free-fall time (see Section 12,21 for details). The equation of 
state (EOS) is shown in column 2, with 'iso' denoting isothermal, 
and 'lar' denoting a piece-wise barotropic EOS (see Section[2]and 
Equations [T] and [2J that is similar to that suggested by Larson 
(2005). The seed column shows which of the two random real- 
isations of the initial turbulent velocity field was used for that 
calculation. N* , M, total an d M* are the total number and mass 
in sink particles and mean sink particle mass, respectively, with 
M, = M t< total / N* . Finally, /m>im q denotes the fraction of 
the total mass in sinks which is made up of sinks with masses of 
1 or greater. 
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Figure 3. The mass in sinks as a function of time for the clouds 
with o^in = ^kin / l^grav | = 0.1). The long-dashed and short- 
dashed lines show the SFE(t) for two isothermal runs which have 
different initial random realisations of the turbulent velocity field. 
The dotted line traces the SFE(f-) for the simulation with the 
barotropic EOS. For comparison, we plot the SFE(t) from the 
isothermal ct^m = 1 simulation (solid line), which shares the same 
seed as the cr^n = 0.1 denoted by the long-dashed line. 



form in the unbound clouds are able to use up their gas reser- 
voir without as much influence from their neighbours. The 
unbound clouds hence move towards a more isolated mode of 
star formation, whereby individual regions of star formation 
evolve independently. 

The second effect which shapes the mass function is the 
degree of fragmentation. From Figure[4] one can see that the 
more unbound simulations still form fairly high-mass sinks. 
What makes the mass function progressively flatter is the 
lack of lower-mass objects, since the clouds with higher lev- 
els of turbulent energy form systematically fewer bound ob- 
jects. This occurs because many regions take several sound 
crossing times to become bound, during which time the in- 
ternal energy has erased much of the structure. Naturally, 
this is more pronounced in the barotropic EOS calculations, 
since the gas heats up as it becomes less dense. The lack of 
fragmentation within the star forming regions reduces the 
role of competitive accretion. 

A s imilar effect is seen in the driven turbulent simula- 
tions of iKlessenl (|200ll ). When the turbulent driving scale 
is small, the stars tend to form in relative isolation. Since 
they have little gravitational interaction with their siblings, 
the resulting mass functions are flat. In contrast, large scale 
driving leads to a highly clustered environment, in which 
competitive accretion is important. Consequently, these lat- 
ter simulations produce mass functions which are more sim- 
ilar to the IMF. 

After considering the unbound clouds, one may expect 
that the clouds with a^in = 10 to result in mass functions 
which are very steep. Figure shows that this is not the 
case. In these clouds, where the collapse proceeds in a quasi- 
homologous fashion, the sinks interact strongly with their 



neighbours, which leads to a mass spectrum that is broadly 
simil a r to the IMF. Indeed , it has been shown l)Zinneckerl 
1 19821 ; iBonnell et all l2001bh that the mass function should 
never get infinitely steep under these conditions, with an 
asymptotic limit for the slope reaching —2.5. Again we see 
that the isothermal and barotropic EOSs give very similar 
results. 

Lastly, we note that when competitive accretion occurs 
in clouds with cv^m ^ 1, the mass function of the accreting 
objects always has roughly the same shape. This in demon- 
strated in Figure [6] which shows the evolution of the mass 
function for one of the isothermal runs with akin — 1- This 
is an important feature of the mass function in these cloud, 
since it means that the star formation need not be termi- 
nated at some special time, but can indeed be terminated 
at any point in the evolution and will still give the correct 
IMF.' 

One can see why the shape of the mass function remains 
constant in the CYkin ^ 1 calculations by considering the 
time-scales for accretio n and fragmen t ation i n these calcula- 
tions. As discussed by IBonnell et all (|2001al ). the accretion 
time-scale in such clouds is comparable to the free-fall time. 
Further, the crossing time, which dictates the decay ti me- 
scale of the turbulence in the gas l|Mac Low et al.lll998l ). is 
related to the free-fall time in these clouds by, 

tff - ( 3n 2r - 1 a Va (4) 

icr \32Gp) iw VlO kin ' 

where r is the radius of the cloud, w rms is the root-mean- 
square of the turbulent velocities, and p is the mean density 
in the cloud. Naturally, this is also the average time-scale 
on which regions of gas will lose their turbulent support and 
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start to collapse, and as such can be thought of as the frag- 
mentation time. Since the time-scales for the accretion and 
fragmentation are similar throughout the cloud's evolution, 
the rate at which new sinks form is then similar to rate at 
which mass is added to the pre-existing sinks via accretion. 
Thus, both the Salpeter section of the mass function and the 
flatter section at lower masses grow together, which main- 
tains the shape of the mass distribution. 

In the contrasting picture to competitive accretion, in 
which the stellar masses are set purely by the fragmenta- 
tion of the molecular cloud, it is n ot clear that the emerging 
IMF w ould be time independent. IClark. Klessen fc Bonnelll 
(|2007l ) pointed out that any stellar mass function that re- 
sults from the collapse of a bound distribution of cores will 
change with time, provided that the time-scale for the cores 
to collapse is a fun ction of t h eir m ass. While this may not be 
a problem for the ISalpeterj (|l955l ) section of the core mass 
function, the time-scale issue may change the mapping be- 
tween the ob served cores and fi nal stars at lower masses (for 
example, see lAndre et al.ll2007f l. 



5 DISCUSSION 

The results of this paper show that arbitrarily low star for- 
mation efficiencies are possible with just slight variations 
in the initial kinematic state of the cloud or star forming 
region. As determinations of the total cloud masses and ki- 
netic energy content are uncertain, such variations are easily 
within our current understanding of molecular clouds. For 
example, an initial o^in = 0.3 corresponds to a mass three 
times less than what is required to be bound or equivalently 
an internal velocity dispersion 1.8 times higher. 

In order to fully understand the significance of 
these results, we need to put them in a broader con- 
text of cloud formation and evolution. Recent work 
has shown that a dynamical formation mechanism 
for molecular clouds is likely due to some larger 
scale compres s ion jlBallesteros-Paredes et al.l 1 19991; 



Hartmann et all l200ll; iBonnell et al.l 120061 ; iDobbs et al l 



20061 ; Vazquez-Semadeni et al. 20061 ) . Shocks induced by 
spiral arms, or supernovae are two such possibilities. In 
these cases, there is no reason to assume that the clouds are 
gravitationally bound on the largest scales. Instead, as local 
regions become cl ose to being gravitat i onally bound due to 
the c ompression l|Klessen et al.l 120051 ; | B all csteros-Paredel] 
120061 ). then star formation is initiated on this local scale 
(Bonnell et al 2006). The star formation is likely to be 
inefficient except in those regions which are strongly 
gravitationally dominated, resulting in a low overall star 
formation efficiency. 

Although the above scenario is dynamical in nature, it 
can easily be mistaken as being s tow in terms of the sta r 
formation rate per free-fall time (|Krumholz fc Tu"]|2007l ). 
The primary reason for this is that an observer is unlikely 
to see only the final stages of the process but is equally 
likely to observe any intermediate state in the evolution. 
Thus what we need to compute is the time-averaged star 
formation efficiency for a given cloud. For example, given 
the end results presented in Figure [5] the time averaged 
star formation efficiency per free-fall time (as defined from 
the initial density) for the cases of akin = 1, 0.5, 0.3 are 



SFE/t ff w 0.18, 0.07, and 0.02, respectively. These are sig- 
nificantly lower than just the final star formation efficiency 
and could be taken to imply that star formation in these sim- 
ulations was slow. Furthermore, we need also consider the 
formation timesca le of the cloud which effectivel y doubles 
the cloud lifetime (|Ballesteros-Paredes et aT1ll999h and thus 
halves the above estimates of the star formation efficiency 
per free-fall time to 0.09, 0.035 and 0.01. 

The results of this paper may also provide an explana- 
tion for the relative lack of star formation in the Madellena 
GMC. We show here that star formation becomes less effi- 
cient and more d istributed as th e clou ds become more un- 
bound. Similarly, IWilliams et al.l (|l994T ) have demonstrated 
that the structure in the Maddelena cloud is significantly 
unbound, and such structures have estimated masses in the 
range from around 2 to 1000 Mq. We hence suggest that the 
highly unbound and dynamic environment in the Madde- 
lena cloud may be responsible for the very low level of star 
formation. 

We also predict that dynamically unbound regions 
should produce more isolated stars that follow a flatter IMF 
than those in the bound clustered environments. Given the 
unbound state of the maj ority of the struct ure in, for exam- 
ple, the Orion A cloud (|Ballv et all I1987T ) or the Rosette 
molecular cloud (|Williams et al.l 19941 ). we would predict 
that a low level of star formation should permeate most 
molecular clouds, providing that the dense regions have a 
sufficient number of Jeans masses. Current observations sug- 
gest that this is indeed the case. The 2MASS and Spitzer 
surveys have revealed distributed popula tions of stars in 
a number of near by star-forming regions (|Carpenterll2000l ; 
lAUen et al.l l2007l) . and the data suggest a smooth tran- 
sition between 'di stributed' and 'clustered' star formation 
l|Allen et al.l l2007t ). From the wide variety in dynamical 
states for the sub-structure in molecular clouds, our results 
would also suggest a smooth transition between heavily clus- 
tered environments and regions of relatively isolated star 
formation. 

One potential problem with the model we present here 
is that the resulting field-star IMF in the Galactic disc would 
become contaminated by the flatter IMFs from the unbound 
regions. If most star formation occurs in the more massive 
bound clusters, then competitive accretion is able to produce 
the observed IMF. In contrast, if the majority of stars in- 
stead form in loose associations, in which gravity only dom- 
inates at the scale of individual stars or systems, then our 
results would suggest an integrated IMF that is shallower 
than the observed distribution. However, the star formation 
efficiency is also important. Although the more unbound 
clouds produce flatter IMFs, they do so with a lower effi- 
ciency, and thus the overall effect of the distributed popu- 
lation is not clear. Estimates for the distributed population 
in several nearby star forming regions suggest that these ob- 
jects comprise around 20 to 30 p ercent of the host cloud's 
star formation (I Allen et al.ll20071 l. which will also include 
any stars that have left dense clusters or groups. This im- 
plies that the contamination proble m is not too severe. We 
also note that lLada fc Ladal (|2003l ) estimate that the ma- 
jority of stars do form in clusters. 

Lastly, we issue a word of caution when interpreting the 
mass functions presented in this study. Sink particles only 
store the potential mass reservoir for the protostar, or pro- 
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Figure 4. The mass functions of the sink particles for the simulations with cski n = i?kin / I -Egrav 5^ 1, given after a time t = 2%. 
The simulations with an isothermal EOS are shown in the upper panels, while those with the barotropic EOS arc shown in the bottom 
panels. The long-dashed distributions show the results from a different realisation of the initial turbulent seed. The isothermal «kin = 0.1 
distribution (filled region) is shown alongside the a^in = 0-3 results. The dotted distribution in each panel is the functional form of the 
IMF for the local field stars given by Chabrier (2003). Note that the mass resolution in all cases is 0.05 Mq, and those sinks with masses 
below this value are still accreting their neighbours. 



tostellar system, and do not allow material to escape once it 
has been accreted. In reality stars will have winds and out- 
flows, which are expected to lower the SFE of the protostel- 
lar cor e. However for low-mass st ars, this effect is likely to be 
small l|Wuchterl fc Klesserfooil ). The formation of higher- 
mass stars may be less efficient due to th e increased radiation 
pressure and presence of ionising winds (|Yorke fc Sonnhalterl 
l2002l ; lKrumholz, McKee fc Kleinll2005l ). Therefore, if the in- 
ternal SFE of a sink particle is a decreasing function of 
increasing mass, then the final stellar IMF would indeed 
be shallower than those presented in this study. To recover 
the observed Salpeter slope, the sink particle mass function 
would then need to be somewhat steeper. The sinks used 
in this study also have acc retion radii that are larger than 
typical binary separations (|Duquennov fc Mayor! [l99ll ). As 
such, they are more representative of multiple systems, for 
systems smaller than the accretion radius r acc (as we point 
out in section [5]). One would therefore expect the mass func- 
tion to steepen somewhat for higher resolution. 

6 SUMMARY 

This paper examines the properties of the star formation 
in a series of 1000 Mq clouds, in which the initial level of 
turbulent support is altered, with akin = Skin / |-Egrav| vary- 
ing from 0.1 to 10. We perform these calculations for two 



different EOSs, first assuming isothermality and then as- 
suming a barotropic EOS, the latter of which is similar to 
that suggested by Larson (2005). The starting density in the 
barotropic EOS clouds is over an order of magnitude lower 
than that in the isothermal clouds. 

We find that a wide range of star formation efficien- 
cies is possible. In the simulations for which a^n ^ 1, the 
star formation efficiencies range from up to 60 percent, to as 
low as around 0.3 percent, after a time of 2 initial free-fall 
times. Since these calculations do not contain any model for 
the feedback from young stars, these SFEs are upper lim- 
its. The isothermal and barotropic equations of state yield 
similar results. We suggest that the very low level of star 
formation in the Maddelena cloud may be a result of the 
highly unbound struc t ure th at has been observed in this 
object IWilliams et~aH (|l994h . In contrast, the simulations 
with akin = 0.1 are able to reach a star formation efficiency 
of around 35 percent, in only 1 free-fall time. 

The mass functions of the 'sink' particles formed in the 
simulations are also found to be related to each cloud's ini- 
tial dynamical state. In the clouds which start with akin ^ 1, 
the mass functions of the sinks becomes progressively flat- 
ter with increasing initial energy in the turbulence. For the 
clouds with akin = 1 to 0.1, the sink particle mass function 
are broadly similar to the stellar IMF. However we caution 
the reader that the sink particles used in this study are larger 
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Figure 5. The sink mass functions for the weakly supported sys- 
tems (those with = i?kin / I -Egrav | = 0.1), given after a time 
t = iff. The solid and dashed distributions show the mass func- 
tion for the isothermal runs, for two different seeds, and the dotted 
distribution is for the simulation with the barotropic EOS. The 
dotted distribution in each panel is the funct i onal form of the 
IMF for the local field stars given bv lChabrieil ll2003h . Note that 
the mass resolution in all cases is 0.05 Mq, and those sinks with 
masses below this value are still accreting their neighbours. 



than the typical binary system, and so are best interpreted 
as a system mass, which would be expected to be flatter 
than a true stellar mass function. All the features seen in 
the mass functions in this study are parallelled in both the 
EOSs that are employed. 

We also point out that for a^in ^ 1, the general shape 
of the mass function is constant in time, including position 
of the turnover. This is an important feature of the mass 
function: regardless of when competitive accretion is inter- 
rupted by gas expulsion, the resulting mass function will be 
consistent with the observed distribution. As such, it does 
no t suffer from the 't ime-scale' problem, discussed in detail 
bv lClarketal] ll2(K)7h. which can arise when trying to map 
a pre-stellar core mass function to a star or system mass 
function. 

We conclude that if competitive accretion is to play an 
important role in shaping the stellar IMF, then the major- 
ity of star formation needs to occur in bound, collapsing 
regions, at the scale of the cluster formation. However, the 
GMCs that host these cluster forming regions need not be 
globally bound entities. Within the GMCs unbound sub- 
structure, we predict that there can be a distributed mode 
of star formation, provided the gas in these structures has a 
sufficient number of Jeans masses. These stars are expected 
to exist in relative isolation, with a trend towards progres- 
sively more isolation in more unbound or dynamic regions. 
We predict that the IMF of this distributed population is 
flatter than that found in either the clustered environment 
or the local field star population. 
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Figure 6. The mass function of sink particles, at four different 
instances, for isothermal simulation with velocity seed number 2 
and akin = 1. The times for which the mass functions arc shown 
are 0.75 (long-dashed), 1 (dot-dashed), 1.5 (short-dashed) and 2 
(solid) initial free-fall times. Note that while the mass function 
at 2 iff is the same as plotted in Figure [4] (top left panel, dashed 
distribution), the binning here is slightly larger, simply to make 
the plot easier to read. Again, the dotted line denotes the Chabrier 
(2003) IMF. 
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